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Abstract 

A specific linear combination of the total solar neutrino fluxes must equal the 
measured solar photon luminosity if nuclear fusion reactions among light ele- 
ments are responsible for solar energy generation. This luminosity constraint, 
previously used in a limited form in testing the no neutrino oscillation hypoth- 
esis, is derived in a generality that includes all of the relevant solar neutrino 
fluxes and which is suitable for analyzing the results of many different solar 
neutrino experiments. With or without allowing for neutrino oscillations, the 
generalized luminosity constraint can be used in future analyses of solar neu- 
trino data. Accurate numerical values for the linear coefficients are provided. 
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I. INTRODUCTION 



If nuclear fusion reactions among light elements are responsible for the solar luminosity, 
then a specific linear combination of solar neutrino fluxes must equal the solar constant. 
More explicitly, we shall see that 

where L is the solar luminosity measured at the earth's surface and 1 A.U. is the average 
earth-sun distance. The coefficient «j is the amount of energy provided to the star by nuclear 
fusion reactions associated with each of the important solar neutrino fluxes, <3?j. Equation ([!]) 
is known as the luminosity constraint. The utility and power of the luminosity constraint 
are due to the fact that the numerical values of the ctj's are determined by the differences 
between nuclear masses and are independent of details of the solar model to an accuracy of 
one part in 10 4 or better. 

The conservation laws of baryon number, lepton number, and charge are the basis for 
the luminosity constraint and are contained within the reaction equation that describes the 
fusion of hydrogen to make helium within stars, 

Ap — > 4 He + e + + 2z/ e + (thermal energy). (2) 

The amount of thermal energy that is delivered to the star depends upon the particular set 
of reactions that occur in realizing Eq. (fj). If, for example, the relation shown in Eq. @ 
involves the production of the relatively high energy 8 B neutrinos, then less thermal energy 
is communicated to the star than if only low energy pp neutrinos are produced. 

The sum in Eq. ([!]) should be taken over all the neutrino fluxes whose associated nuclear 
fusion reactions could in principle contribute significantly to the energy budget of the sun. 
If no detailed knowledge derived from solar models is used to limit the appearance of terms 
in Eq. ([[]), then the luminosity constraint can be applied to tests of the hypothesis of no 
neutrino oscillations. The luminosity constraint provides an additional condition that must 
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be satisfied by neutrino fluxes that are otherwise allowed to have arbitrary amplitudes when 
fit to the available solar neutrino data PHIO]. 



A. Previous work 

Spiro and Vignaud ]I| first proposed, in a lucid and insightful paper, the use of the 
luminosity constraint as a test of the null hypothesis for solar neutrino propagation, i.e. as 
a test independent of solar models of the assumption of no neutrino oscillations. Following 
these authors, most of the pioneering applications of the luminosity constraint (see, e.g., 
@-^|,||) have approximated the solar neutrino spectrum by grouping the neutrinos into three 
sets, the low-energy (principally pp )neutrinos, intermediate energy neutrinos (usually taken 
to be either 7 Be neutrinos only or the 7 Be neutrinos plus the CNO and pep neutrinos), and 
finally the high-energy ( 8 B) neutrinos. This approximation was necessary and appropriate 
when the number of solar neutrino experiments was small (two or four), but is no longer 
required or optimal now that the number of solar neutrino experiments is six (chlorine ||11||, 



Kamiokande 0, SAGE 0, GALLEX + GNO 130, Super-Kamiokande 0, and SNO 
fl7p and growing (BOREXINO H] and ICARUS PI), with additional experiments in the 
planning stages. 

The derivation given here explains (and, in some cases, corrects) the the results that 
were simply stated in Ref. |§. 

All previous discussions with which I am familiar, including Ref. |J, have implemented 
the luminosity constraint in the context of showing that solar fluxes with arbitrary ampli- 
tudes, but without the distortion of the energy spectrum implied by neutrino oscillations, do 
not fit the available data on the measured rates of solar neutrino experiments. The papers 
by Hata, Bludman, and Langacker [||], Parke ||, Heeger and Robertson ||, and Bahcall, 
Krastev, and Smirnov ]T(]] were important in persuading many physicists who are not fa- 
miliar with solar models that a particle physics solution was required for the solar neutrino 
problem. 
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B. What is this paper about? 



My goal in this paper is to provide a general formulation for the luminosity constraint 
that can be used in future analyses, with or without allowing for neutrino oscillations, that 
may include six or more experiments. I also want to provide a specific derivation, lacking in 
the literature, for the coefficients in Eq. ([!]). The lack of a general derivation in the literature 
has led to a confusion about the basis for the luminosity constraint and to significant errors 
in the published values of some coefficients. 

In this paper, I derive the coefficients for the luminosity constraint for all seven of the 
important neutrino fluxes shown in Table U After deriving the coefficients for the general 
form of the luminosity relation, I discuss the most appropriate approximations to make in 
analyzing data sets in which the number of measured neutrino event rates is not sufficient to 
allow a statistically meaningful application of the full luminosity constraint. The discussion 
in this paper is a natural generalization of the treatment given in the very important paper 
by Hata et al. 0, which presented cogently the argument that the measured solar neutrino 
event rates (chlorine' ||11|| , Kamiokande ||12|| , SAGE [0], and GALLEX [|14|]) required new 



physics, even before the epochal Super-Kamiokande [16| and SNO |Ty|] measurements. 

The formulations presented in previous discussions can be recovered from the present 
analysis by assuming that one or all of the following assumptions is valid: i) certain neutrino 
fluxes are zero; ii) the CNO neutrino fluxes (from 13 N and 15 beta-decay) are equal; or iii) 
the standard solar model ratio of pep neutrino flux to pp neutrino flux is correct. 

In the future, the generalized luminosity constraint can and should be implemented in 
analyses that determine solar neutrino parameters. The additional constraint provided by 
the measured solar luminosity will be especially important when pp and 7 Be neutrino fluxes 
are measured as well as the 8 B neutrino flux. As more experimental data become available, 
the analyses of neutrino oscillations will become more independent of the standard solar 
model and it will be natural and convenient to incorporate the luminosity constraint, Eq. (|l]). 

The luminosity constraint can be written conveniently in a dimensionless form by con- 
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sidering the ratios of all the neutrino fluxes to the values predicted by the standard solar 
model, and dividing both sides of Eq. ([I]) by Lq/[4tt(A.U.) 2 ). We obtain 

1 = ? ( Toiev) ' (3) 
where the dimensionless neutrino fluxes, (pi, are the ratios of the true neutrino fluxes to the 



neutrino fluxes predicted by the BP2000 standard solar model pOj, i.e. 



<f>i = $,/$ i (BP2000) . (4) 

The quantities are the ratios of $j(BP2000) to the characteristic solar photon flux defined 
by L /[47r(A[/.) 2 (lOMev)]. Thus 

Oi = $i(BP2000)/ (8.5272 x 10 10 cm" 2 s" 1 ) . (5) 

In calculating the characteristic solar flux, I have used the recent best-estimate solar lu- 
minosity (see Ref. |21|), 3.842 x 10 33 ergs^ 1 , that is derived from all the available satellite 
data. 

Depending on the context, I shall use to refer to either the neutrino flux produced 
locally or integrated over the entire sun. This dual usage will not cause any confusion since 
the specific meaning of (pi will be clear in all cases. Moreover, because of the linearity of the 
averaging process, equations that are valid locally have the same form when the results are 
integrated over the entire sun. 

Table | provides the numerical values for the dimensionless form of the luminosity con- 
straint. 

Some people may prefer to use the luminosity constraint in its dimensional form (Eq. p. 
In this case, one can ignore the values of aj given in Table | and use just the tabulated values 
of a,. If one wants to use the constraint in its dimensionless form but with a different basis 
set of neutrino fluxes instead of those from BP2000, then one can replace the aj given here 
by new values computed from Eq. (|5|) with the different model fluxes. 
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TABLE I. Luminosity constraint: neutrino characteristics. The average neutrino ener- 



gies are taken from Ref. |23| for 7 Be and from Ref. [24| for all other sources. The neutrino energies 
include thermal effects from electron and ion motion and from the solar temperature profile, as 
well as atomic ionization effects. The nuclear data are taken from p3]. The quantities a and a are 
defined in Eqs. (H)-(H), the dimensionless form of the luminosity constraint. 



r nix 


Reaction 


v^w© 










(MeV) 


(10 MeV) 




4>(pp) 


p + p — > 2 H + e + + v e 


0.2668 


1.30987 


0.6978 


(f>(pep) 


p + e~ + p -> 2 H + v e 


1.445 


1.19193 


0.001642 


<p(hep) 


3 He + p -> 4 He + e + v e 


9.628 


0.37370 


1.09E-07 


0( 7 Be) 


7 Be + e~ -> 7 Li + 


0.814 a 


1.26008 


0.05594 


<K 8 B) 


8 B -> 8 Be + e+ + i/ e 


6.735 


0.66305 


0.0000592 


0( 13 N) 


13 N _^ 13 C + e + + „ e 


0.706 


0.34577 


0.006426 


0( 15 O) 


15 -» 15 N + e+ + i/ e 


0.996 


2.15706 


0.005629 



a 89.7% 0.8631 MeV and 10.3% 0.3855 MeV. 

The linear equality of Eq. (|3|) must be supplemented by the physical requirement that the 
number of nuclear reactions that terminate the proton-proton chain not exceed the number 
of initiating nuclear reactions, which can be translated into the following inequality [f|, 

( 7 Be) + ( 8 B) < 0(pp) + 0(pep) . (6) 

The physical basis of Eq. (^j is that the 3 He nuclei, which ultimately give rise to 7 Be and 8 B 
neutrinos via the nuclear reaction 3 He(a,7) 7 Be, are created by pp and pep reactions. One 
pp or pep reaction must occur in order to supply the 3 He nucleus that is burned each time a 
7 Be or 8 B neutrino is produced. In principle, Eq. (§) (or Eq. HD considered separately permits 
a 'Be neutrino flux that is twice as large as is allowed by Eq. (H). Since the 14 N(p, 7) 15 
reaction is the slowest process in the CNO cycle, one must also have 

0( 15 O) < 0( 13 N) . (7) 

With the currently available solar neutrino data, the additional constraints provided by 
Eq. (||) and Eq. ([7]) are automatically satisfied by the best-fit solutions for the undistorted 
neutrino fluxes to the measured event rates p2| . 

In Sec. [TIJ I derive explicit expressions for the a coefficients that appear in Eq. (^) in terms 
of measured atomic mass differences and computed neutrino energy losses. I summarize in 
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Sec. |TTT| the principal assumptions that are used in the derivation and discuss the results and 
application strategies in Sec. [TV]. 



II. DERIVATION OF THE LUMINOSITY CONSTRAINT 

In this section, I shall derive expressions for the energy coefficients, c^, that appear in the 
luminosity constraint, Eq. ([!]). Section [11 A| treats the simpler case of the CNO neutrinos and 



Sec. |11 ij| derives the coefficients for neutrino fluxes produced by reactions in the pp chain. 
The numerical values given in Table | were calculated using Eqs. (|TU]) and ( |TTD of Sec. [LI A 
and Eqs. ©-PD of Sec. ITTE. 



The coefficients aj that appear in the luminosity constraint are independent of details 
of the solar model to high accuracy. As we shall see below, the a's represent the differ- 
ences between nuclear masses with only small corrections, included here, due to the average 
thermal energy of the fusing particles. The largest of the thermal corrections amounts to 
0.7% of the average neutrino energy loss for the pp neutrinos pi] and 0.15% for the 7 Be |23 



neutrinos. For all other solar neutrino sources, the thermal corrections to the energy loss are 
much smaller |26[] . Since the total neutrino energy loss is itself only a small fraction of the 
a's we compute in the following two subsections, one can show that thermal energy affects 
the values of a given in Table | by less than or of order 0.01% for all of the important solar 
neutrino energy sources. 

In carrying out the calculations, we will use Rij to represent the reaction rate per unit 
of time per unit of volume between two fusing ions, i and j, where 

(i, j)n(i)n(j) 

Here is the local thermal average of av , the product of the relative velocity of particles 
i, j and their interaction cross section, and the Kronecher delta prevents double counting of 
identical particles. For example, -R34 = ( 3 He, 4 He)n ( 3 He) n ( 4 He) and R pp = (p,p)n(p) 2 /2. 

In the following, I shall denote the average neutrino energy from a particular nuclear 
reaction, X, by < E v > (X). Table | lists accurate values for these neutrino energies. The 



values given in Table | are averaged over the energy spectrum from each source and include 
corrections for solar effects such as contributions from the thermal motion of the fusing ions, 
averages over ionization states, and the temperature profile of the sun. The average solar 
neutrino energy losses are taken from Refs. The masses, for example M ( 13 C) or 

M ( 1 H), that appear in the equations of Sec. [11 A] and Sec. [11 B| are atomic masses which I 



have taken from Ref. 25 



For convenience in the calculations, we will introduce a fictitious neutrino flux density, 



$33 = ( 3 He, 3 He)n( 3 He, 3 He)/2 



(9) 



which is the rate of the 3 He( 3 He, 2p) 4 He reaction (and would be the flux density of neutrinos 
produced by this reaction if the 3 He — 3 He reaction gave rise to neutrinos). This fictitious 
flux will not appear in any of the final formulae. 



A. CNO neutrinos 

The 13 N beta-decay corresponds to the thermal energy derived from the reactions 
12 C(p, 7) 13 N and 13 N -> 13 C + e+ + v e . Hence 

a ( 13 N) = M ( 12 C) + M ( 1 H) - M ( 13 C) - (E„) ( 13 N) . (10) 

Similarly, the energy associated with the 15 O beta-decay derives from the reactions 
13 CQo,7) 14 N, 14 N(p,7) 15 0, 15 -> 15 N + e+ +u e , and 15 N(p,a) 12 C. Therefore, 

a ( 15 0) = 3M ( 1 H) + M ( 13 C) - M ( 4 He) - M ( 12 C) - (E v ) ( 15 0) . (11) 

The neutrino flux from 17 F beta-decay is a potential measure of the primordial 16 
abundance in the sun (see Ref. |27j ) , but does not play a significant role in the generation of 
the solar luminosity. For completeness, I include here the coefficients a and a that describe 
the reactions 16 0(p, 7) 17 F and 17 F — > 17 + e + + v e . The appropriate linear coefficients for 
use in Eq. (^]) are 



S 



a 



( 17 F) = 2.363 MeV , a ( 17 f) = 1.09 E - 07 . (12) 



The value of a ( 17 F) is so small (because of the high Coulomb barrier for this reaction) that 
the flux 0( 17 F) is not relevant for practical applications of the luminosity constraint. 

B. pp neutrinos 

The analysis of reactions in the pp chain is simplified by the fact that 2 H and 3 He are 
burned quickly at solar temperatures [ . The lifetime for nuclear burning of 2 H is ~ 10 -8 



yr and the lifetime of 3 He is ~ 10 5 yr. These values are both very small compared to the 
10 10 yr lifetime of a proton (which is destroyed primarily by the pp reaction). Therefore, it is 
an excellent approximation to assume that both 2 H and 3 He are in local kinetic equilibrium 
(rate of destruction equals rate of production). 

From the local equilibrium of 2 H, dn( 2 H)/dt = 0, one has 

( 1 H, 1 H)rt (^rl) 2 /2 + ( 1 H,e- + H)n(e)n(H) = ^H, 2 H)n ( x H) n ( 2 H) . (13) 



Equation fll3|) states that the production of deuterium via the pp and pep reactions is 
balanced by the destruction of deuterium via the 2 H(p, 7) 3 He reaction. The equilibrium of 
3 He, dn( 3 B.e)/dt = 0, implies 

(% 2 R)n ( 1 H) n ( 2 H) = ( 3 He, 3 He)n(3) 2 + ( 3 He, 4 He)n ( 3 He) n ( 4 He) 

+ ( 3 He, 1 H)n ( 3 He) n ( x h) . (14) 



Equation ([14]) describes the fact that the rate of production of 3 He via the 2 H + 1 H reaction 
is balanced by the destruction of 3 He via the 3 He — 3 He, 3 He- 4 He, and hep reactions. For the 
term describing the 3 He- 3 He reaction in Eq. (|14"D, the factor of one-half from the identity of 
the fusing particles is cancelled by the factor of two representing the destruction of two 3 He 
ions. 

Combining Eqs.(|), and (0), we find 

$(p-p) + $(pep) = 2 x $(3, 3) + <f>{hep) + $ ( 7 Be) + $ ( 8 b) . (15) 



Let £j represent the thermal energy released to the star as a result of the nuclear fusion 
reactions associated directly with each neutrino producing reaction. Then we can write 

L Q / Atc(A.U.) 2 = CNO terms + e pp ®(p-p) + e pep $(pep) 

+ (e 33 /2) [$(p-p) + ®(pep) - <$>(hep) - $ ( 7 Be) - $ ( 8 B)" 

+ e hep $(hep) + (e 34 + e e7 ) $ ( 7 Be) + (e 34 + e p>7 ) $ ( 8 B) . (16) 

I have used Eq. (|15|) to eliminate the fictitious flux 0(3,3) from Eq. (|16|) . This substitu- 
tion associates with each of the real neutrino fluxes in the pp chain an additional energy 
contribution proportional to e 33 , the energy released to the star via the 3 He- 3 He reaction. 
Physically, these terms proportional to e 33 represent a way of keeping track of how much 
energy from the the 3 He- 3 He reaction should be associated with the other neutrino fluxes. 

The values of e« can be calculated by writing explicitly the reaction equations that are 
associated with the production of each neutrino flux. One finds 

e pp = 3M ( 1 H) - M ( 3 He) - (E u )(p-p) , (17) 

e pep = 3M ( 1 H) - M ( 3 He) - (E v )(pep) , (18) 

e 33 = 2M ( 3 He) - M ( 4 He) - 2M ( 1 h) , (19) 

e 34 = M ( 3 He) + M ( 4 He) - M ( 7 Be) , (20) 

e e7 = M ( 7 Be) + M ( 1 h) - 2M ( 4 He) - (E v ) ( 7 Be) , (21) 

e P:7 = M ( 7 Be) + M ( 1 H) - 2M ( 4 He) - (E v ) ( 8 B) , (22) 



and 



e hep = M ( 3 He) + M ( 1 h) - M ( 4 He) - (E v )(hep) . (23) 
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In calculating e e 7, one must average over the two 7 Be neutrino lines with the appropriate 
weighting and include the 7-ray energy from the 10.3% of the decays that go to the first 
excited state of 7 Li. 

The values of the a's can be determined using the following relations between the a 
coefficients and the e coefficients that follow from Eq. ([H]). We have 

a(p-p) = e pp + 0.5e 33 , a (pep) = e pep + 0.5e 33 , a ( 7 Be) = e 34 + e e7 - 0.5e 33 , 

a ( 8 B) = e 34 + e pJ - 0.5e 33 , a(hep) = e hep - 0.5e 33 . (24) 



III. WHAT ASSUMPTIONS ARE MADE IN DERIVING THE LUMINOSITY 

CONSTRAINT? 

The basic assumption made in deriving the luminosity constraint is that nuclear fu- 
sion reactions among light elements are responsible for the observed solar luminosity. More 
specifically, I assume in Sec. |I| that the specific nuclear reactions that have been recog- 
nized [f28H 32|l over the six decades since Hans Bethe's epochal work on the subject as being 
most important at temperatures of order a keV are indeed the fusion reactions that power 
the sun. The characteristic temperature of the sun can be estimated relatively well without 



making use of a detailed model [28,29 



In order to derive Eq. ( |T3"D and Eq. (|H"D , it is necessary to assume that both 2 H and 3 He 
are in local kinetic equilibrium (rate of creation equal rate of destruction). As discussed 
in Sec. |11B| , this is an excellent approximation because the lifetimes for nuclear burning of 
these isotopes are short compared to the evolutionary time scale for the sun. 



IV. DISCUSSION 

I have given in Sec. [IT] an explicit derivation of the luminosity constraint and have pre- 
sented in Table | coefficients that can be used in the dimensionless form of the constraint, 
Eq. (HD . The coefficients that are given in Table | are calculated accurately and include small 
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corrections to the neutrino energy release that result from the high temperatures (~ keV) in 
the region in which fusion reactions occur in the sun. The thermal corrections to the energy 
release are less than or of order 0.01% in all cases. 

The form of the luminosity constraint given here includes all the important solar neutrino 
fluxes. One can recover approximately the coefficients used by previous authors who have 
combined neutrino fluxes, or who have considered only a reduced set of fluxes, by making 
the relevant choices among the fluxes listed in Table |. However, the reader is warned not 
to expect precise agreement; there are many inaccurate numerical values in the published 
papers.Q 

The generalized form of the luminosity constraint presented here can be used, as the more 
restricted constraint has been used in the past, to help test the validity of the null hypothesis 
for solar neutrino oscillations. In this test (cf. Refs. (l|-0), the neutrino fluxes are allowed to 
have arbitrary amplitudes subject to Eq. (|3]) and the condition that the energy spectra are 
undistorted by neutrino oscillations. One should also impose the two inequalities, Eq. (|6|) and 
Eq. (0), that follow from the requirement that the number of nuclear fusion reactions that 
terminate the proton-proton chain and the CNO cycle not exceed the number of inititiating 
nuclear reactions. 

In the past, applications of the luminosity constraint have been limited to tests of the no 
oscillation hypothesis. Nothing in the derivation of the luminosity constraint (or the supple- 
mentary inequalities, Eq. ^ and Eq. [7]), requires this limitation in the range of applications. 
In the future, when more experimental data are available, the luminosity constraint can be 



x In fact, I am responsible for an egregious but unimportant error. In Ref. Q, the listed value for 
a{hep) is wrong because I neglected to subtract the value of €33/2 = 6.4298 MeV in going from e^ep 



to a(hep) [see Eq. (24) of the present paper]. The error is unimportant, although embarrassing, 
because the Super-Kamiokande experiment has placed fl6| a strong upper limit on the hep flux; 



this upper limit is also in agreement with the predicted hep flux for the standard solar model [20], 
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used together with the measured solar neutrino interaction rates, energy spectra, and time 
dependences, to help determine neutrino oscillation parameters. 
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